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Abstract. In this paper we consider an equilibrium last-passage percolation model on an envi- 
ronment given by a compound two-dimensional Poisson process. We prove an L^-formula relating 
the initial measure with the last-passage percolation time. This formula turns out to be a useful 
tool to analyze the fluctuations of the last-passage times along non-characteristic directions. 



1. Introduction and the main result 

1.1. The last-passage percolation model. Let P C be a two-dimensional Poisson random 
set of intensity one. On each point p G P we put a random positive weight ujp and we assume 
that {cjp : p E P} is a cohection of i.i.d. random variables, distributed according to a distribution 
function F, which are also independent of P. Throughout this paper we wih make the fohowing 
assumption on the distribution function F of the weights: 



(1) / e"^ dF{x) < +00 , for some a > . 

^0 



'0 

This condition was used in [7] to prove the existence of invariant measures for the Hammersley's 
interacting fluid process we will introduce below. For each p, q G M^, with p < q (inequality in 
each coordinate, p / q), let n(p, q) denote the set of all increasing (or up-right) paths, consisting 
of points in P, from p to q, where we exclude all points that share (at least) one coordinate with 
p. So we consider the points in the rectangle ]p, q], where we leave out the south and the west side 
of the rectangle. The last-passage time between p < q is defined by 



L(p,q) := max { V Wp'} 



Ti7en(p,q) 

p ero 

When F is the Dirac distribution concentrated on 1 (each point has weight 1 and we will denote 
this F by 5i), then we refer to this model as the classical Hammersley model [ll[9j. 

A crucial result is the following shape theorem (see Theorem 1.1 in [8], p. 164): set = (0,0), 
n = {n, n), 

(2) 7 = 7(F) = sup ^^^^"'"^^ >0 and f{x,t):=-fV^. 
Then ^{F) < oo and for all x, t > 0, 

(3) hm ^IMlZl!:!)) = li^ WMrx^ ^ ^(^^ _ 

r— )-oo T r—^oo f 

1.2. The interacting fluid system formulation. It is well known that the classical Hammersley 
model has a representation as an interacting particle system [U |9] . The general model has a similar 
description, although a better name might be an interacting fluid system. We start by restricting 
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the compound Poisson process {up : p G P} to M x M_|_. To each measure z/ on M we associate a 
non-decreasing process z/(-) defined by 

( ] - S ^((O'S;]) for X > 
- I -i/((x,0]) for X < 0. 

Let J\f be the set of all positive, locally finite measures v such that 

liminf ^ > 0. 

y^-oo y 

We need this condition to define the evolution of the process, since otherwise all mass will be 
pulled to minus infinity. The Hammersley interacting fiuid system [M^ : t > 0) will be defined as 
a Markov process with values in AA, as was done in [7j. Its evolution is defined as follows: if there is 
a Poisson point with weight a; at a point {xo,t), then M^{{xq}) = M^_{{xq}) and for x > xq, 

Maixo,x]) = {MUixo,x])-oo)+. 

Here, Mi_ is the "mass distribution" of the fiuid at time t if the Poisson point at {xQ,t) would be 
removed. To the left of xq the measure does not change. In words, the Poisson point at {xQ,t) 
moves a total mass u to the left, to the point xq, taking the mass from the first available fiuid to 
the right of xq. (See Figured] for a visualization, in case of atomic measures, of the process inside 
a space-time box.) 

In this paper we follow the Aldous and Diaconis [1] graphical representation in the last-passage 
model (compare to the result in the classical case, found in their paper): For each u £ M, x G M 
and t > let 

(4) L^{x, t) := sup {u{z) + L{{z, 0), {x, t))} . 

The measure defined by 

M^{{x, y]) := L^{y, t) - L^{x, t) for x < y , 
defines a Markov process on M and it evolves according to the Hammersley interacting fiuid system 

We now make the following important observation for a random initial condition u, which basi- 
cally follows from translation invariance. 

Theorem 1. Suppose v £ M is a random initial measure on M independent of the Poisson process 
in R X M+, whose distribution is translation invariant. For any speed V G M and any x G M, we 
have 

Ly{Vt,t) = L^{x,t) - u{x - Vt). 

The relevance of this result is most clear when we consider equilibrium measures of the Ham- 
mersley 's interacting fiuid process. Assume that we have a probability measure defined on M and 
consider u G M as a realization of this probability measure. We say that v is time invariant for the 
Hammersley interacting fiuid process (in law) if 

Mi' = M^ = v for all i > . 

In this case, we also say that the underlying probability measure on N is an equilibrium measure. 
It is known that there is only one family of ergodic equilibrium measures for the Hammersley 
interacting fiuid system [7]. Let us denote it by {i^a : > 0}) where 

(5) \:=¥.ux{l). 

For simple notation, put L\ := Ly^. The main result of this paper is the following formula: 
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Figure 1. In this picture, restricted to [0, x], the measure consists of three atoms 
of weight 5, 3 and 7. The Poisson process, restricted to [0,x] x [0,t], has two points 
with weights 4 and 7. The measure Mj^^^ consists of three atoms of weight 1, 4 and 
6, while at time t, it consists of one atom with weight 7. A total weight of 4 + 6 has 
left the box due to Poisson points to the left of the box, while a total weight of 2 
has entered. 



Corollary 1. Recall ([2]) and ([5]), and let 

(6) "^-©^ "-"i'^'-'^y 

Here, Vx is the characteristic speed corresponding to Lx and tpx is the growth rate of Lx{Vxt,t). 
Then 

(7) E( {Lx{x, t) - [ux{x - Vxt) + V'At]}' ) = Var (La {Vxt, t) ) . 

1.3. A central limit theorem for the classical model. To illustrate the importance of ([7]), let 
us restrict ourselves to the classical Hammersley model. In this set-up, the equilibrium measures 
are one-dimensional Poisson processes of intensity A, and 7 = 7(5i) = 2. Thus, 

1 2 
Vx ■= j2 A ■ 
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Cator and Groeneboom [6] proved that the variance of Lx grows sub-hnearly along the characteristic 
speed A^^. Together with Corollary [H this implies 

Corollary 2. Let {zt)t>o be a deterministic path. Then 

(8) 1-^ mL,iz,,t)-[i.,izt-X-H) + 2X-H]f) _ ^.^ yar{L,{X-\t)) 

t— >oo t t— i>oo t 

Proof of Corollary [21 Formula ([7]), applied to the classical model, gives us 

E( {Lx{zt, t) - [uxix - X-H) + 2X-h] }' ) = Var(LA {X-\ t) ) . 
On the other hand, ^ shows that 

,.^Var(MA-M))^„_ 

t— >CJO t 

which proves (j8]). □ 

Corollary 3. Let {zt)t>o be a deterministic path such that 

lim — = a . 

t—>oo t 

Then 

(9) lim ^ = a := \aX — —\ . 

t^oo t X 

Furthermore, if a ^ A^^ then 

(10) lim ¥(Lx{zt, t)<Xzt + l + {aVt)u) = ¥{N < u) , 
where N is a standard Gaussian random variable. 

Proof of Corollary [3} Corollary [2] shows that 

Lx{zt,t) - [ux{zt - X-H) + 2X-h] 
lim = , 

i-i-oo 

in the sense. Since vx is a one-dimensional Poisson process of intensity A, this implies Q and 
(HOD. □ 



Remark 1. Cator and Groeneboom [6] proved that yYav(^Lx (A ^t,t) ) is of order t^^^ , which gives 
us the same order for the iJ'-distance between Lx{zt,t) and [z^a(3; — A^^t) + 2A^^t] . 

Remark 2. The central limit theorem for Lx (along any direction) was proved by Baik and Rains 
[5j . Their method was based on very particular combinatorial properties of the classical model that 
do not seem to hold for the general set-up. Our approach reveals the strong relationship with the 
initial configuration. 

Remark 3. In the general set-up, Corollary{l\ implies: If the variance of Lx along the characteristic 
speed Vx is sub-linear, and the equilibrium measure has Gaussian fluctuations, then Lx will also 
have Gaussian fluctuations along non- characteristic directions. 
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Remark 4. For the classical Hammersley process an important formula for the variance of Lx{x, t) 
was derived in [6j, Theorem 2.1: 

Var(LA(x, t)) = -Ax + ^ + 2AE(x - Xx{t))+, 

A 

where X\{t) is the position at time t of a second class particle starting at zero. This formula was 
pivotal in deriving the cube-root behavior of Lx in [6j, and later corresponding formulas were used to 
prove cube-root behavior for TASEP [3] and for ASEP However, this formula does not directly 
show the relationship with the initial configuration. Also, there seems to be no direct way to deduce 
([7|) from this formula, even if we reformulate it, as was done in Equation (3.6) of [6\, in terms 
of the exit-point of the longest path from (0,0) to {x,t), which is the right-most z for which the 
supremum in @ is attained. 

2. Proof of Theorem [T] and Corollary [T] 

Recall that 

L^{x, t) = sup Mz) + L{{z, 0), (x, t))} . 

z<x 

Clearly, 0), (T^t, t)) = L{{z + x — T^i,0), {x,t)). By assumption, v has a translation invariant 

distribution, independent of L. This implies that {z i— t- v{z)} = {z i— t- v{z + x — Vt) — v{x — l^t)}, 
and 

Li,{Vt,t) = suY> {v{z + x-Vt)-u{x-Vt)+L{{z + x-Vt,d),{x,t))} 

z<Vt 

= sup {v{z) + L{{z, 0), (x, t))} - v{x - Vt) 

z<x 

= Li,{x,t) -u{x - Vt). 

This proves Theorem [TJ □ 

Corollary [1] now follows from results in [7]: there it is shown that for any speed V , the stationarity 
of Lx leads to 

¥.Lx{Vt,t) = VXt + ^-i'^t/X. 

This follows from the fact that the Hammersley fluid process has intensity A on the bottom side 
of the rectangle between (0,0) and (x,t), and intensity 7^/(4A) on the left side (this refers to the 
expected mass of the fluid leaving the interval [0, x] through per time unit). When we define the 
characteristic speed Vx = 7^/(4A^), then 

¥.Lx{Vxt,t) = i^xt. 

This together with Theorem [T] immediately shows ([7]). □ 



3. The lattice last-passage percolation model 

In the lattice last-passage percolation model one considers i.i.d. weights {wp : p £ 7?}, dis- 
tributed according to a distribution function F. For F{x) = 1 — e~^ (exponential weights), we have 
a similar shape theorem as ^ with limit shape given by 

f{x,t) = {^+Vt)^. 
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We know from |3j that the invariant measures are given by 

y 

z=x+l 

where {X^ : z G Z} is a collection of i.i.d. exponential random variables with parameter p. The 
analog to formula ([7]) is 

E( {Lp{x, t) - [vp{x - [Vpt\ ) + i,pt]Y ) = Var (Lp (L^ptJ , t) ) , 

where 

p2 



Vp ■= J- r^- and 



Together with the cube-root asymptotics this implies that 

E( {Lp{zt,t) - [upizt - Vpt) + V]}' ; 



Therefore, if 



then 



and if a ^ Vp then 



lim 

t—>-oo t 



lim — 

t— >oo t 



lim = a — 



lim P [Lpizt^t) <- + + {aVt)u \ = ¥iN < u) 

■^oo y p I — p J 



t— >oo 

where is a standard Gaussian random variable 



Remark 5. Ferrari and Pontes [10] determined the dependence on the initial condition for the 
totally asymmetric exclusion process, which is isomorphic to the lattice last-passage percolation 
model with exponential weights. The method developed in this paper resembles the ideas in their 
paper. Baldzs [2j used a different method to get a generalization of the Ferrari- Pontes result for 
certain types of deposition models. It is not clear to us whether our methods would work for these 
more general deposition models. 

Remark 6. In the general lattice model, the shape theorem ([3]) holds. However, not much is known 
about the limit shape f . If this function would not be strictly curved (we know it is convex, so this 
would mean that there are "flat" pieces), then the methods used in ^ to prove the existence and 
uniqueness of semi-infinite geodesies in a fixed direction do not apply, and we are not able to prove 
the existence of equilibrium measures. 
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